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ABSTRACT 

A long infrared jet has been discovered by the Spitzer c2d legacy program 
in core A of L1251. It is associated with a very embedded Class object with 
an accretion luminosity of about 0.9 L & derived by radiative transfer model 
fitting to the observed SED. Comparing the observed IRAC colors along the 
infrared jet with those calculated from a model of an admixture of gas with a 
power-law temperature distribution indicates that the jet is possibly created by 
a paraboloidal bow shock propagating into the ambient medium of ^(H 2 ) = 10 5 
cm -3 . In addition, the variation of the power- law index along the jet suggests 
that the portion of hot gas decreases with distance from the jet engine. The 
molecular outflow in this region has been mapped for the first time using CO 
data. From the calculated outflow momentum flux, a very strong lower limit to 
the average accretion luminosity is 3.6 L & , indicative of a decrease in the 
accretion rate with time. 

Subject headings: ISM: individual (L1251) — stars: formation 
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1. INTRODUCTION 

Stars form out of the collapse of dense molecular gas cores. However, this process is 
complex, and the details are far from understood. Young stellar objects (YSOs) and their 
immediate surroundings are expected to evolve significantly throughout star formation. For 
example, YSOs will evolve in luminosity and temperature, as their masses grow and they 
evolve toward the main sequence (e.g., Myers et al. 1998). In addition, the associated gas 
will evolve dynamically. As collapse proceeds gas flows from cores through disks and onto 
protostars, driving jets and outflows in the process. These outflows are a crucial mechanism 
by which angular momentum is removed from a star forming system where accretion occurs 
through a rotating disk. Therefore, an outflow is a good tracer of star formation, especially 
when the central protostar is deeply embedded. 

Large molecular outflows have been mapped in (sub) millimeter lines of CO (Fukui 1989, 
Morgan et al. 1991), which trace the matter cooled down after being shocked. However, more 
recently shocked material such as warm molecular hydrogen produces ro-vibrational or pure 
rotational transitions at infrared wavelengths. The Infrared Spectrograph (IRS) aboard the 
Spitzer Space Telescope (SST) has been used to map the shocked H 2 line emission (Neufeld 
et al. 2006, Maret et al. 2009). Lately, Neufeld & Yuan (2008) showed that the Spitzer 
Infrared Array Camera (IRAC) intensities were dominated by warm molecular hydrogen 
emission in the supernova remnant IC 443. In addition, Ybarra & Lada (2009) also showed 
that IRAC colors could be used to study properties of the hot, shocked gas components in 
HH54. 

The Spitzer Legacy Program "From Molecular Cores to Planet Forming Disks" (c2d; 
Evans et al. 2003) discovered a new jet feature on the western edge of L1251 (core A in 
Sato et al. 1994), whose distance is 300(±50) pc (Kun & Prusti 1993). No outflows were 
previously detected in this region although there are three IRAS sources (Kun & Prusti 1993; 
Sato et al. 1994). No optical jets or Herbig Haro objects have been reported in this region, 
probably due to the high optical depth. The discovered infrared jet is the longest (~ 1 pc) 
except HH 111, which is onesided and ~ 8 pc (Reipurth 1989, Reipurth et al. 1997). In this 
paper, we examine the properties of warm molecular gas along the jet feature with the IRAC 
images. We also study the molecular outflow and the embedded YSO associated with the 
jet using, respectively, CO data and a set of infrared through millimeter photometric data. 
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2. OBSERVATIONS 

The Spitzer c2d Legacy program observed core A of L1251 (L1251-A) at 3.6, 4.5, 5.8, 
and 8.0 /j,m with the Infrared Array Camera (IRAC; Fazio et al. 2004) on 4 September 
2004 and at 24 and 70 fim with the Multiband Imaging Photometer for Spitzer (MIPS; 
Rieke et al. 2004) on 30 November 2004 (PID:139; AOR keys: IRAC 5165824 and MIPS 
9425664). The detailed description of the data processing can be found in the c2d final 
delivery documentation (Evans et al. 2007). 

To cover the whole region of the infrared jet, the CO 2—1 transition at 230.537970 
GHz has been mapped with the 6 m telescope at the Seoul National Radio Astronomical 
Observatory in March and April 2009. The beam FWHM is 48" at 230 GHz. The main beam 
efficiency and pointing accuracy are 0.57 and ~ 3", respectively. The velocity resolution is 
0.127 km s -1 after binning by two channels. 

For the analysis of the Spectral Energy Distributions (SEDs), we also use 350 fim and 
1.2 mm continuum emission maps observed with the Submillimeter High Angular Resolution 
Camera II (SHARC-II) mounted on the Caltech Submillimeter Observatory, and the Max- 
Planck Millimeter Bolometer (MAMBO) at the IRAM 30m telescope, respectively. Details 
of these continuum data can be found in Wu et al. (2007) and Kauffmann et al. (2008). 

3. RESULTS 

Figure la shows a three-color image comprised of IRAC 3.6, 4.5, and 8.0 /im images. 
The dominant feature in the image is a jet extending about 10'. The extended blue color 
emission perpendicular to the jet is consistent with the distribution of the millimeter dust 
continuum emission (see upper panel of Figure lb), so it is likely dominated by scattered 
light within a dense cloud with some contribution from shocked gas associated with jets and 
outflows. There is also a small bipolar nebula structure, similar to that seen in core E of 
L1251 (L1251B, Lee et al. 2006), to the east of the long jet. 

In Figure la, stellar objects (IRS3 and IRS4) associated with the jet and bipolar nebula 
are not clearly seen. However, these central engines of the jet and bipolar structure are well 
detected at 24 and 70 /im (see Table 1 for fluxes). At 70 /im, the two infrared sources become 
brighter (see upper panel of Figure lb). In addition, as seen in Figure lb, these sources are 
each located at both 350 /im and 1.2 mm dust continuum peaks, indicating they are likely 
very embedded objects. In L1251-A the c2d team identified 4 YSOs, marked as IRS1, IRS2, 
IRS3, and IRS4 in Figure lb. The infrared source between IRS3 and IRS4 is classified as a 
galaxy based on the c2d criteria (Evans et al. 2007). Table 1 lists fluxes of the four YSOs 
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at various wavelengths, and the SED of the YSO associated with the jet (IRS3) is shown in 
Figure 2 (as red diamonds). The SED clearly indicates that IRS3 has a very thick envelope 
with strong (sub)millimeter emission. IRS4 is associated with the bipolar nebula and shows 
a similar SED to that of IRS3, although it is not presented in this paper. 

Figure lb (lower panel) shows the CO 2—1 outflow map on top of the IRAC 4.5 /im 
image. The outflow map covers the whole jet feature as well as the bipolar nebula. Along the 
jet, it shows two distinct lobes with little overlap, suggesting a large inclination of rotational 
axis with respective to the line of sight. The outflow lobes are well correlated spatially with 
the infrared jet. No clear outflow structure is seen towards IRS4 and its associated bipolar 
nebula, although some weak red-shifted emission is present. 

4. ANALYSIS 

The jet feature associated with IRS3 appears consistently in all IRAC bands, indicating 
that the jet emission is probably produced by a common gas component such as molecular 
hydrogen, which has its line transitions distributed over all IRAC bands. Therefore, we 
calculate the IRAC colors of [5.6 /xm]/[8.0 /zm] and [4.5 /xm]/[8.0 /im] to compare with those 
of thermal gas modeled to explain the molecular hydrogen emission from the jet (Figure 3). 
The colors were calculated in 4 pixel (~ 5") radius apertures at the positions marked with 
black circles in Figure 2. In the model calculation of the colors (Figure 3), the molecular 
gas was assumed to be composed of H2 and He, with n(He) = 0.2 71(^/2) ■ We calculated 
the level population of H 2 gas in statistical equilibrium. The collisional coefficients were 
obtained from Le Bourlot et al. (1999), and an ortho-to-para ratio of 3.0 is assumed. The 
temperature of the shocked H 2 gas was modeled in two ways: isothermal and power-law 
admixed. In the latter, the infinitesimal H2 column density has a power-law relation with 
the temperature T, dN ~ T~ b dT, in the range of 300 ~ 4000 K. More detailed descriptions 
can be found in Neufeld & Yuan (2008) and Shinn et al. (2009a, b). 

We know that this jet is located inside a molecular cloud (Sato et al. 1994), thus, 
the extinction close to the jet engine is possibly much higher than 0.5 mag, which is the 
interstellar extinction toward L1251-A (Kun & Prusti 1993). Therefore, the expected colors 
are extinction corrected with Ay ~ 5 mag, which is obtained from the 1.2 mm dust continuum 
emission at 0.1 pc from IRS3, following the "Milky Way, R v = 3.1" curve (Weingartner & 
Draine 2001, Draine 2003). Extinction corrected colors with Ay ~ 0.5 and 5 mag are not 
different. However, close to YSOs (C and E in Figure la), Ay ~ 50 mag. We tested this 
high extinction only for the colors of C and E, which move diagonally up and right in the 
color-color diagram. 
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According to the comparison between the observed and modeled colors along the jet, 
the IRAC colors can be explained with the power-law admixture model, as Neufeld & Yuan 
(2008) and Shinn et al. (2009a, b) found. All colors are located between the model of 
n(iJ 2 ) — 10 6 cm -3 and the LTE line, where the pre-shock density must be higher than 
10 5 cm -3 (Wilgenbus et al. 2000). There might be some contribution from other emission 
components not considered in our model, such as the CO v — 1 — transitions in IRAC 2. 
However, such contamination is not likely according to Neufeld & Yuan (2008). The power- 
law index, b is the smallest (~ 2) at the jet feature close to the central embedded source 
(C in Figure la) and the bipolar nebula feature (E in Figure la) in the east of the long jet, 
indicating that there are more hot components close to the embedded YSOs. Alternatively, 
the scattered light in the outflow cavity may enhance the colors (more scattered light at 
shorter wavelengths). If we apply Ay = 50 mag, the colors at C and E become similar 
to those of the model with b~ 1. Except for these two points, b is between 3.3 — 5.2, 
consistent with results of Neufeld & Yuan (2008), who showed that this range of b values 
can be explained by the paraboloidal bow shock model. Bow shocks can be considered 
as a geometrical summation of planar C-shocks, whose postshock temperature is nearly 
isothermal (Neufeld et al. 2006). In this sense, Neufeld & Yuan (2008) showed that a single 
bow shock generates a shocked H 2 gas with b~ 3.6, and b can be higher than 3.8 when slow 
bow shocks, where H 2 is not dissociated even at the apex, are mixed with fast bow shocks. 

Following Young et al. (2004), Bourke et al. (2006), and Dunham et al. (2006), we 
modeled the SED of the YSO associated with the jet (IRS3) using the 1-D dust continuum 
radiative transfer code DUSTY (Ivezic et al. 1999). In addition to a central stellar object, a 
1-D disk model (Butner et al. 1994) was adopted with a surface density profile of S(r) oc r~ L5 
and temperature profile of T(r) oc r~ a35 to calculate the internal input SED for the best-fit 
model. The envelope is assumed to have a power-law density profile and is heated by the 
interstellar radiation field attenuated by Ay =3 mag. For the dust opacity, we adopt the 
same dust model used by Crapsi et al. (2008): a mixture of carbon (29% of the total dust 
mass) and silicate grains covered with various molecular ices. For a detailed explanation of 
the modeling procedure see Young & Evans (2005) and Dunham et al. (2006). The best-fit 
SED to the observed fluxes is presented in Figure 2. The best-fit model has an internal 
luminosity of ~ 0.9 L Q with T star = 2500 K. This internal luminosity derived by the model, 
which is considered to be the true bolometric luminosity of the source, is greater than 
calculated from observed fluxes because some of these long-wavelength fluxes were measured 
with apertures smaller than the total extent of the source. The envelope of the model has a 
power-law density profile index of p = 1.55, a mass of ~ 12 M , and inner and outer radii 
of 250 and 18,000 AU, respectively. The internal luminosity derived by the model can be 
considered the current accretion luminosity. 
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The average accretion luminosity over the YSO lifetime can be derived from the momen- 
tum flux in the outflow. We use the CO J= 2 — 1 map to calculate the outflow momentum 
flux, Fqq, following the method presented by Hatchell et al. (2007). This quantity is the av- 
erage rate at which momentum is injected into the outflow, and we calculate Fqo > 5.9 x 10~ 6 

M km s _1 yr _1 and list this result in Table 2. There is an inclination dependence that 
we leave in the result since the source inclination is not known. Two additional unknowns 
in the above calculation are the optical depth and the temperature of the outflowing gas 
(needed to convert line intensity to column density assuming LTE). Following M. Dunham 
et al. (2009, in prep), we assume the gas is optically thin and has a temperature of 17.6 
K, which minimizes the line intensity to column density conversion. Optically thick gas and 
any other temperature in the range of 10 — 100 K will increase Fqo, thus the above result is 
a very strong lower limit to the true value. 

Since molecular outflows are driven by the transfer of momentum from a jet/wind ejected 
by the protostellar system to the ambient medium and the jet/wind ejection process is closely 
tied to accretion onto the protostar, the calculated Fqo, which represents the average rate 
at which momentum is injected into the outflow, is related to (M acc ), the time-averaged 
protostellar mass accretion rate (see Bontemps et al. 1996): 



Jcnt M w V w 

In the above equation, M w is the mass-loss rate in the jet/wind, V w is the velocity of the 
jet/wind, and / ent is the entrainment efficiency between the jet/wind and the ambient gas. 
With typical values for the above parameters of = 0.1, V w ~ 150 km s^ 1 , and / ent = 0.25 

J" acc 

(Bontemps et al. 1996 and references therein), we calculate (M acc ) > 1.6 x lO -6 ^^ M q 
yr _1 . An outflow dynamical time (r^ > 14.2 x 10 4 ^| yr) is calculated by dividing the extent 
of the outflow (400"/sini) by the outflow velocity (~ 4/cosi km s _1 ). Assuming accretion at 
the above rate for this length of time gives an accreted protostellar mass of M acc > 0.22 
M . Finally, accretion at the above rate onto a protostar with this mass and R = 3 R gives 
a spherical accretion luminosity of L acc = GM ^ ac ^ > 3.6 £ . 



5. DISCUSSION 

We calculated the bolometric luminosity and temperature of each YSO identified in 
L1251-A, following the method used in Dunham et al. (2008) (see Table 1). The apparent 
drivers (IRS3 and IRS4) of the infrared jet and bipolar nebula have < 1 L Q , similar to 
IRS2 in L1251B (L1251B-IRS2), another source associated with infrared bipolar structure 
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(Lee et al. 2006). Although IRS3 and IRS4 are ~ 10 times fainter than L1251B-IRS2 at 24 
/iin, they are very bright at wavelengths longer than 70 /jm. Their bolometric temperatures 
are ~20 K, much lower than that of L1251B-IRS2 (140 K). Finally, the ratio of L smm /Lf, i 
for IRS3 and IRS4 (~0.13 and 0.17, respectively, where L smm is the luminosity at A > 350 
/iin) is much greater than that of L1251B-IRS2 (~0.03), indicative of Class sources with 
very thick envelopes. 

The dominant feature in L1251-A is the long, well-collimated jet associated with IRS3. 
This type of jet structure can be produced in the IRAC bands by ro-vibrational or pure 
rotational line emission of H2 in a pulsating jet with little precession (model Al in Smith 
& Rosen 2005). According to Smith & Rosen (2005), this jet model (Al) transfers the bulk 
kinetic energy to large distances but excites H 2 with the lowest efficiency However, the jet 
shows undulating structure toward both ends. Smith & Rosen (2005) showed this kind of 
helical stream in a jet with a slow uniform-speed precession. Therefore, the jet in L1251-A 
might be explained with a small precession angle (but greater than that of Al) and a long 
pulsating period (much longer than 60 years). In fact, the jet associated with the northern, 
red component of the CO outflow features two distinct emission bulks, indicative of episodic 
ejection. The time interval between episodes of ejection is 1300/sini years if we assume a 
jet speed of 150 km s -1 . 

The infrared jet and the CO outflow are well correlated in shape and length as seen in 
Figure lb. From the CO 2 — 1 map, we calculated a time-averaged mass accretion rate and the 
accretion luminosity expected from accretion at this rate (Table 2). In spite of uncertainties 
in / ent and M w /M acc , the calculated time-averaged accretion rate (and resulting accretion 
luminosity) are strong lower limits, as explained in §4. In addition, a large inclination (there 
is little overlap between red and blue components of the CO outflow and no significant 
intensity difference between two infrared jet lobes) of 70° will give L acc ~ 83 L , which is 
much greater than L acc derived from dust modeling (0.9 L Q ). Since the luminosity from dust 
modeling depends on the current mass accretion rate, this rate must be lower now than the 
time-averaged value over the outflow lifetime. 

The comparison between the observed IRAC colors and the ones from the model of a 
power-law admixture of gas (see Figure 3) indicates the amount of hot components varies (in 
general, more hot gas closer to the jet engine). If we use the same colors ([3.6]-[4.5] vs. [4.5]- 
[5.8]) as those presented by Ybarra & Lada (2009), who modeled only the hot (2000—4000 
K) gas component, C and E (b ~ 2, see Figure 3) are located at n(H) > 10 5 cm" 3 , while 
N2, N3, and S (b ~ 3.3 — 4.5) lie close to the model line of n(H) = 10 5 cm" 3 . Therefore 
the density inferred by their model is not very different from that inferred from our model, 
although, unlike in our model, Ybarra & Lada (2009) included collisions between H 2 and H. 
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However, Nl is placed at n(H) ~ 2 x 10 3 cm -3 in their model, consistent with our model 
result since the model of Ybarra & Lada (2009) considered only the hot gas component, and 
in our model, b ~ 5.2 for Nl suggests a relatively small amount of hot gas at that position. 
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Table 1. Properties of YSO candidates in L1251-A 



FLUXES (mjy) 
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DEC. 
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Lbol (Lq) 
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IRS4 f 
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75 13 37.1 


0.43 


1.07 


0.75 


0.37 


1.90 
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7500 
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0.6 


21 


0.13 



a (SHARCII) aperture size = 
b (MAMBO) aperture size = 
c SSTc2d J222933.4+751316 
d SSTc2d J222959.5+751403 
e SSTc2d J223031.9+751409 
f SSTc2d J223105.6+751337 
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80 " 
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Table 2. Outflow and Mass Accretion Properties 



Quantity Units Value 

Outflow Momentum Flux (Fco) 
Average Protostcllar Mass Accretion Rate ((M acc )) 
Dynamical Time (r d ) 
Protostcllar Mass Accreted (M acc ) 
Accretion Luminosity (L acc ) 



a i is the inclination of the rotational axis with respect to the line of sight. 



KT b M km s- 1 yr- 1 > 5.9 
1CT 6 Mq yr- 1 > 1.6 

10 4 yr > 14.2 C %1 

■ 

— cos % 



Mq > 0.22 

L« > 3.6^2- 




Fig. 1. — (a) Three-color composite Spitzer image of L1251-A. The IRAC 3.6, 4.5, and 
8.0 fim data are presented, respectively, as blue, green, and red. Green lines are placed 
along the jet and bipolar nebula features. Black circles show the regions where the colors 
plotted in Figure 3 were derived, and white circles to the east and west of each black circle 
indicate where the background was measured, (b) Upper Panel: 1.2 mm MAMBO intensity 
map (contours) on top of MIPS 70 /zm image. IRS3 and IRS4 are located at the peaks of 
MAMBO emission. Another emission peak exists between IRS1 and IRS2, which might be 
a prestellar core. Lower Panel: CO molecular outflow map (contours) on top of IRAC 4.5 
/im image. The dotted line encloses the total area mapped, and the green filled circle at the 
right bottom denotes the beam at 230 GHz. The blue contours are integrated from -13.0 
to -7.0 km s _1 , while the red contours are integrated from -1.0 to 5.0 km s~ x . The contours 
start at 1.2 and 1.6 K km s _1 for blue and red components, respectively, and increase by 0.4 
K km s _1 . The lobes of the main outflow are well correlated with the infrared jet along the 
NS direction. 
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Fig. 2. — The best-fit model SED (green solid line) to the observed fluxes (red diamonds) of 
IRS3, with the internal input SED shown (blue dashed line). Blue diamonds represent model 
fluxes in the apertures used for photometry. Crosses indicates the SED emerging only from 
the envelope. This SED modeling gives an internal luminosity for IRS3, produced mostly by 
accretion, of ~ 0.9 L . 
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Fig. 3. — Two-color diagram of [4.5 /xm]/[8.0 jum] and [5.6 /xm]/[8.0 /xm] expected from pure 
H2 emission with an assumption of A v = 5 mag. Blue lines: single-component models for 
n(H 2 )= 10 4 (dot-dot-dot-dashed), 10 5 (dot-dashed), 10 6 (dashed), 10 7 cm~ 3 (dotted), and 
LTE (solid). Open squares represent the results of a range of temperatures up to 2000 K 
(top right of each curve) in steps of 100 K. Red lines: models with a power-law temperature 
distribution, with dN = aT~ b dT and b in the range of 2 (top right of each curve) to 6. 
The line type for each model is the same as in the single- component models. Each symbol 
represent the color measured at each position marked in Figure la. The error bars include 
the flux measurement error and the flux calibration error (~ 10%). 



